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The notion of a non-perturbative effect is ambiguous if it requires the subtraction of a perturbative 
part defined by a diverging series. A common procedure consists in dropping the order of mini- 
mal contribution and the higher orders. This allows us to isolate very accurately the one-instanton 
effect for the double-well potential. For the one plaquette gauge theory, an exact analytical expres- 
sion can be written for the non-perturbative part. We report recent attempts to extend this approach 
to the average plaquette of quenched QCD. Our goal is to express the non-perturbative effects in 
terms of expressions of the form (f5) B e~ A P calculable semi-classically. The situation is compli- 
cated by zeroes of the partition function in the complex j3 plane (presumably near 5.75 ± i0.2). 
We discuss two methods to describe the intermediate and large order behavior of the perturbative 
series. One is inspired by mean field theory (logarithmic specific heat) and reproduces accurately 
the known perturbative series with only two free parameters. A diagrammatic interpretation of 
this fact is still lacking. The other is based on infra-red renormalons with a factorial growth show- 
ing up at order larger than 20 and a possible effective theory interpretation. These extrapolations 
are compatible with the non-perturbative part of the plaquette being proportional to a 4 . We pro- 
pose an exponential parametrization to the corrections to the universal part of the beta function 
and find results compatible with the suggestion of a 2 corrections made by C. Allton. 
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In these proceedings, we will discuss the definition of the nonperturbative part of the average 
plaquette and its parametrization in terms of expressions of the form A/3 B e~ c ^. There is a reason- 
able case that this is feasible for Wilson action in the fundamental and the hope is that A, B and C 
could be calculable semi-classically. Technical details regarding recent progress are available in a 
recent preprint [jl|]. Instead of trying to cover most of the material discussed in this preprint, we 
will rather emphasize selected points that we have disccused with some of the participants during 
the conference. 

We would like to emphasize that the commonly used "rule of thumb" for perturbative series (at 
a given coupling, you drop the order in perturbation theory that gives the smallest contribution and 
all the higher orders), is a poor man substitute for regularizing the diverging perturbative series by 
introducing a large field cutoff or a large action cutoff as recently proposed [§, ||]. The field cutoff 
can then be fixed by an optimization procedure based on the strong coupling expansion [Q, HJ. 
However, a first step consists in understanding the error associated with this simple procedure. For 
a generic asymptotic series A ~ ^k a kp k , we can define the error at order k: 

k 

Afe(jS) =A nwnerica i(P) -Y,aip-' . (1) 

/=o 

If we assume that ~ f}~ k ~ l ak+i , (for /3 large enough) and that the large order behavior is 
\ak\ ~ | |C2|^r(A: + C3), we find ^ that the minimal error is 

Min k \A k \ * V^ldKlC^/jS) 1 / 2 -^^ . (2) 

This equation works well for the ground state of the anharmonic oscillator where the exponential 
suppression can be interpreted semi-classically. For the double-well, it only works if we take the av- 
erage of the two lowest energy states in order to cancel the one-instanton effect. The one-instanton 
is a non-perturbative effect larger than the error estimated from the the large order behavior of the 
perturbative series. 

The separation between the perturbative and the non-perturbative part is a computational com- 
modity. There is no claim that this separation has a physical meaning, unless one can find observ- 
ables that are zero in perturbation theory, for example the difference between the two lowest energy 
states for the double-well potential. A simple example [^] where one can see exactly what is in the 
non-perturbative part is the SU (2) one plaquette model 

Z(J8) = / dUc-^- L 2 ReTrU ^=Z Pert XI5)+Z NPert XP) . (3) 

JSU{2) 

Due to the compactness of SU (2), there is no large field problem. It is possible to expand Z as 
a converging sum of /3-dependent coefficients, denoted A*(2j3) below, multiplying fi~ k . We can 
get the regular perturbative series by decompactifying the integrals, namely adding the tails of 
integration from 2/3 to 00, in order to get j8 -independent coefficients. At first sight this amounts to 
neglecting e -2 ' 3 effects, however, this affects the large order of the perturbative coefficients, which 
now grow at a factorial rate. In order to minimize the perturbative error, we truncate the series at 
order r(k*) the integer closest to the one obtained with the procedure discussed above (see Ref. 
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[jl|, ||] for details). More explicitly, 

z Per ,(P) = {M-" 2 ^' 2 X>H/r* , (4) 

Jt=0 

with 

The nonperturbative part consists in the dropped terms minus the added tails: 

Z NPM .(P) = W)-T(P)) t (6) 

oo 

R^) = {^l 2 2 1 ' 2 £ A k (2f5)l3- k , (7) 

r(*»)+l 

m = (^)- 3 /2 2 i/ 2 x ^p-k ^m r dte -^ . ( 8) 

k~o k\[l/2-k)J 2 p 

A detailed calculation shows that at leading order, R and T are both proportional to /3~ 3 / 2 e~ 2 ^, 
while their difference is proportional to /3 2 e 2 ^ as predicted by the asymptotic behavior of the 
regular perturbative series. One can improve the accuracy of the result by keeping a finite range of 
integration for the coefficients and optimizing the range using the strong coupling expansion. If the 
order is large enough, it essentially amounts to use the bound (t < 2/3 ) provided by the compactness 
of the group [Q]. 

The same program can in principle be carried on for quenched QCD in 4 dimensions, however, 
keeping track of multiple tails of integration is a non-trivial task. We believe that the best approach 
is to use the stochastic perturbation method [||] with boundaries or periodic boundary conditions 
in configuration space. Before attempting this calculation for quenched QCD, we have considered 
models where independent numerical calculations are possible, namely one variable integrals and 
ground states of simple quantum mechanical systems [f7[ gjp. In absence of boundary, stochastic 
perturbation theory works quite well for these models, however, introducing a boundary is a delicate 
procedure in part because results may depend on few configurations. M. Naides has investigated 
three methods, reflection at the boundary, rejection of updates going out of the boundary and a- 
posteriori elimination of configurations outside of the boundary and concluded in favor of the third 
method. We are planning to compare these methods for quenched QCD. 

In the following, we will restrict the discussion to the accuracy of the regular perturbative 
series of the average plaquette in quenched QCD (N = 3, D = 4) 

P(fi) = (1/JT P ) h(l - (l/N)ReTr(U p ))\ . (9) 

We used a series in j8 1 calculated ^ up to order 10. A figure with coefficients up to order 16 
is also available Extrapolated estimators [[It]] indicate that P« (1/5.74 — 1/jS) 1 ' 08 in good 
agreement with a ratio analysis [|TT|]. This type of behavior is not expected. The power series in 
j8 _1 should have a zero radius of convergence because the plaquette changes discontinuously at 
j3 — > ±oo [12]. In addition, a singularity exactly on the real axis is not seen in 2d derivative [ |10| ] of 
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Figure 1: Zeroes of the partition function in the complex j3 plane for a 8 4 lattice using reweightings at j3o = 
5.55 (red), 5.6 (cyan), 5.65 (magenta), 5.7 (green) and 5.8 (blue). The dots correspond to distinct bootstraps 
and the solid lines to the successive radius of confidence 



P. Such a singularity would requires massless glueballs which are not known to exist for the pure 
fundamental action. A plausible explanation [[K]] is that the singularity is slightly off the real axis. 
If we assume an approximate logarithmic behavior for the specific heat as in mean field theory 

-ap/ay3«in((i/j3 m -i/y3) 2 +r 2 ), qo) 

we obtain peaks and series compatible with what we know provided that 5.7 < j3„, < 5.9 and 
0.001 < r < 0.01. This suggests the existence of zeroes of the partition function in the complex j8 
plane with 

0.03 - 0.001 fa < |Imj3| < O.OljS^ ~ 0.33 (11) 

A. Denbleyker has performed direct searches for such zeroes using reweigthing methods. On a 
4 4 lattice, he found a pair of zeroes near /3 = 5.55 ± £0.1, clearly within the region of confidence and 
in good agreement with the results on 4L 3 lattices of Ref. [13]. We are not aware of existing results 
on larger isotropic lattices. On a 8 4 , we found no clear evidence for zeroes with |Im/3 1 < 0.03. Our 
preliminary findings are shown in Fig. 1 . The radii of confidence have been drawn without taking 
into account the autocorrelations which should slightly reduce these radii. Unlike the 4 4 case, there 
are no zeroes clearly within the respective radii of confidence. As the radius of confidence increases 
like the square root of the logarithm of the number of independent configurations, it will be difficult 
to find direct evidence for zeroes with |Imj3| « 0.1. Consequently, it would be useful to introduce 
a positive adjoint term in the action in order to move the zero closer to the real axis. If this works, 
we could follow the zero as ^adjoint is decreased and try to extrapolate to fiadjoint = 0. 

A pair of complex singularities near /3 = 5.8 is about all we need to account for most of the 
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known terms of the perturbative series. Integrating the mean field ansatz of Eq. (10), we obtain 



P P ert = £M^-^(Li 2 (i8- 1 /(i8 m 1 +Jr))+h.c, (12) 

k=0 

with U 2 (x) = Li^o-^A 2 • We fixed F = 0.003 and obtained K = 0.0654 and j3 m =5.787 using the 
known values of ag and «io. As shown in Ref. [jl|], this simple model agrees well with the numerical 
values from order 3 to 16. This fact should have a simple Feynman diagram interpretation. 

A pair of complex singularities cannot govern the asymptotic behavior of a series that must 
have a zero radius of convergence. At this point, the best we can do is to use what is known in the 
continuum about the Borel transform of the series Jl4[ [T5| ] and assume that it can be connected to 
the lattice Jl6| ] by a coupling redefinition j3 = j3(l+rfi //$ + ...): 

V = E^« f 2 dte-P'(l-t33/l6n 2 y l - 204 / m . (13) 

k=Q Jt ^ 

Note that ?i = corresponds to the UV cutoff, t% = 16ft /33 to the Landau pole. In the continuum, 
it has been argued [ |T5] ] that it is necessary to introduce the gluon condensate in order to keep ?2 low 
enough and regularize the perturbative series, exactly as we advocated in lattice gauge theory. In 
the following, we take ? 2 = 00 in order to get a regular perturbative series (using the one-plaquette 
analogy). This model has a minimal perturbative error of the form 

Min k \A k \ ~ s.SGS) 204 / 121 - 1 /^ 16 * 2 / 33 ^ . (14) 

Except for the -1/2 in the exponent this is the two loop RG invariant. This is reminiscent of the 
y3 1 / 2 found in the one plaquette case. Following this analogy further, Eq. (14) could in principle 



be compared with what would be obtained from the probability distribution for one plaquette after 
integrating over all the other links. 

The two models defined by Eqs. ( |I2| ) and (13) yield similar coefficients up to order 20-25. 



Above order 25, the renormalon model has larger coefficients growing at a factorial rate. D. Du has 
been trying to put together the two types of behavior in the framework of dispersion relations. These 
two models are compatible [|l|] with the idea already put forward [ 1 1 , ^] that the non-perturbative 



part of the plaquette scales like a 4 with a(/3) defined with the force method §TJ\, |I8[ ]. 

The question is now: can we parametrize a{fi) in a way that would be suggestive from a semi- 
classical point of view? For this purpose, we start with a Taylor expansion [ |l7| ] of a(/$) expressed 
in units of ro = 0.5 f in the interval 5.7 < j3 < 6.92: 

ln(a/r ) = -1.6804 -1.7331 (ft - 6) +0.7849 (fi -6) 2 - 0.4428 (/3 -6) 3 . (15) 

We propose the parametrization ^ 

^ln(a/r )/Jj8 = -(47T 2 /33) + (51/121) j 8" 1 +Ae^ . (16) 



As shown in Fig. 0, in the interval 5.9 < j8 < 6.4, the derivative of Eq. (15) can be fit using 



A = —1.35 10 7 and B = 2.82. It should be noted that when j8 > 6.4, there is a significant difference 



between the derivative of Eq. (15) and the derivative of a similar expression provided in a previous 



paper [ 18]. The difference is hardly visible if we plot ln(a/ro) in the interval 5.6 < j8 < 6.6, but 
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dln(a/r_0) /d/3 c 
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Figure 2: Left: d\n(a / rg) / d j3 using Ref. [ |l7| | (thick dashes), Ref. JT^ ] (small dashes) and our parametriza- 
tion (solid red line). Right: estimation of the (approximate) constant C in Eq. (16). 



when we take the derivative, the large constant term disappears and the cubic term that varies the 
most between the two expressions is amplified by a factor 3. This explains the large differences 
observed in Fig. || when j8 > 6.4. It would be interesting to try to get more accurate Taylor 
expansions in this specific region. 

We can now integrate and fit the constant of integration denoted C: 

ln(a/r ) =C-(47T 2 /33) j 8 + (51/121)ln( j 8)-(A/B)e-^ . (17) 

Fig. |2| shows that C = 4.528(1) provides a good fit in the region where we have good agreement 
with both Refs. JT7| , [l8| ]. A nice plateau appears between j8 = 5.9 and 6.3. The extremal values 
in this interval are 4.5272 and 4.5282. It is also possible to obtain the lattice scale Al from the 
constant of integration C, namely Al = exp(— C)/ro ~ 4.4 MeV. It is possible to use Eq. (|l7|) to 



predict ln(a/rn) at large j3. for instance, at j8 = 7.5, we obtain -3.59 (for -3.63 in Ref. Ql9|]) and 
-4.74 at J8 = 8.5 (for -4.81 in Ref. ^) 

Note that the value of B seems consistent with the idea [ ]2"0| ] of using a 2 pert corrections for 
this quantity. The symbol a pert refers to the one-loop or two-loop expression which in the short 
/3 interval considered here can hardly be distinguished from each other. The assumption of a 2 pert 
corrections fixes B = Sn 2 /33 ~ 2.4 which is close to the value 2.82 obtained above. 

There remain many challenges: check the exponential parametrizations more accurately at 
large j8, calculate the parameters determined numerically semi-classically, find the precise connec- 
tion between the non-perturbative part of the plaquette and what is called the gluon condensate in 
the continuum and find improved perturbative methods with a smaller error in the crossover region. 

We acknowledge valuable discussions at the conference with many participants and in partic- 
ular C. Allton, M. Creutz, F. di Renzo, A. Duncan and A. Velytsky. This research was supported 
in part by the Department of Energy under Contract No. FG02-91ER40664. A. Denbleyker was 
supported by a Van Allen research grant. M. Naides was supported by a REU grant NSF PHY- 
0353509. 
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